1 



(N 
O 
O 
(N 

X 



Modeling Deep Inelastic Cross Sections in the Few GeV Region 
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We present preliminary results on simple modifications and corrections to GRV94 leading order parton distri- 
bution functions such that they can be used to model electron, muon and neutrino deep-inelastic scattering cross 
sections at low energies. (Presented by Arie Bodek at NuIntOl, Dec. 2001, KEK, Tsukuba, Japan). 
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Data from atmospheric neutrino experi- 
ments [|] have been interpreted as evidence for 
y T oscillations with sin 2 2a > 0.88 and 



x 10" 3 < Am 2 < 4 x 10" 



eV 2 . These 
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neutrino data are in the few GeV region. There- 
fore, good modeling of and cross sections at 
low energies is needed. The modeling of and 
cross sections is even more important for the 
more precise next generation neutrino oscillations 
experiments. 

These include MiniBooNE, MINOS, CNGS, 
and experiments in the new neutrino facility to 
be constructed at the JHF high intensity 50 GeV 
proton synchrontron in Japan. 

The quark distributions in the proton and neu- 
tron are parametrized as Parton Distribution 
Functions (PDFs) obtained from global fits to 
various sets of data at very high energies. These 
fits are done within the theory of Quantum Chro- 
modynamics (QCD) in either leading order (LO) 
or next to leading order (NLO). The most impor- 
tant data come from deep- inelastic e//x scattering 
experiments on hydrogen and deuterium, and 
and Vfj, experiments on nuclear targets. In pre- 
vious publications we have compared the 
predictions of the NLO MRSR2 PDFs to deep- 
inelastic e/fi scattering data 10 on hydrogen and 
deuterium from SLAC, BCMS and NMC. In or- 
der to get agreement with the lower energy SLAC 
data (down to Q 2 = 1 GeV/c 2 ), and at the high- 
est values of a; (x = 0.9), we find that the follow- 
ing modifications to the MRSR2 PDFs must be 
included. 



1. The relative normalizations between the 
various data sets and the BCDMS_system- 
atic error shift must be included 



2. Deuteron binding corrections need to be ap- 
plied as discussed in ref. ||. 

3. The ratio of d/u at high x must be increased 
as discussed in ref. || . 

4. Kinematic higher-twist originating from 
target mass effects || are very large and 
must be included. 

5. Dynamical higher- twist corrections are 
smaller but also need to be included [^],[| . 

In addition, our analysis including QCD Next 
to NLO (NNLO) terms shows || that most of 
the dynamical higher-twist corrections that are 
needed to fit the data within a NLO QCD analy- 
sis originate from the missing NNLO higher order 
terms. A recent calculation Q also shows that dy- 
namic higher twist corrections are very small. If 
most of the higher-twist terms that are needed to 
obtain agreement with the low energy data actu- 
ally originate from target mass effects and miss- 
ing NNLO terms, then these terms should be the 
same in and e//x scattering. Therefore, low 
energy data should be described by the PDFs 
which are fit to high energy data and are modified 
to include target mass and higher-twist correc- 
tions to describe low energy e//z scattering data. 
With this idea in mind, we find what modifica- 
tions are need to be applied to GRV94 || leading 
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Figure 1. Electron and muon data on protons 
(SLAC, BCDMS and NMC) compared to the pre- 
dictions with GRV94 PDFs (LO, dashed line) and 
the modified GRV94 PDFs (LO+HT, solid line). 



Figure 2. Electron and muon data on deuterons 
(SLAC, BCDMS and NMC) compared to the pre- 
dictions with GRV94 PDFs (LO, dashed line) and 
the modified GRV94 PDFs (LO+HT, solid line). 



order PDFs such that the PDFs describe both 
high energy and low energy electron/muon data. 

In order to describe low energy data down to 
the photoproduction limit (Q 2 — 0), and account 
for both target mass and higher twist effects, we 
find that the following modifications are required 
for the GRV94 LO PDFs. 

1. We increase the d/u ratio at high x as de- 
scribed in our previous analysis The 
corrections to the u and d distributions are 
described in detail in appendix A. 

2. Instead of the scaling variable x we use the 



scaling variable x w = (Q 2 + B)/(2Mv + A) 
(or =x(Q 2 + B)/(Q 2 + Ax)). This modifica- 
tion was used in early fits to SLAC data |tJ . 
The parameter A provides for an approx- 
imate way to include target mass effects 
and higher twist effects at high x, and the 
parameter B allows the fit to be used all 
the way down to the photoproduction limit 

(Q 2 -o). 

3. In addition as was done in earlier non-QCD 
based fits Q to low energy data, we mul- 
tiply all PDFs by a factor Q 2 / (Q 2 +C). 
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Figure 3. Comparison of representative CCFR 
Vfj, and Vn charged-current differential cross sec- 
tions (flJII on iron at 55 GeV and the predictions 
of the GRV94 PDFs with (LO+HT, solid) and 
without (LO, dashed) our modifications. 



This is done in order for the fits to describe 
both intermediate- x data and data in the 
photoproduction limit, where the structure 
function F2 is related to the photoproduc- 
tion cross section according to 



47t 2 qem 

Q 2 



F 2 = 



0.112to& GeV 2 



Q 2 



4. Finally, we freeze the evolution of the 
GRV94 PDFs at a value of Q 2 = 0.24 (for 
Q 2 < 0.24), because the GRV94 PDFs are 
only valid down to Q 2 = 0.23 GeV/c 2 . 

In our analysis, the measured structure func- 
tions are corrected for the BCDMS systematic 
error shift and for the relative normalizations be- 
tween the SLAC, BCDMS and NMC data @,|. 
The deuterium data are also corrected for nuclear 
binding effects 

A simultaneous fit to both proton and deuteron 
data yields A=1.735, B=0.624 and C=0.188. Fig- 
ures and I show the SLAC, BCDMS and NMC 
data as compared to the predictions of the stan- 
dard leading GRV94 PDFs (LO, dashed line) and 
with our modifications (LO+HT, solid line). 



Our value of C=0.188 is smaller than in other 
analyses || to low Q 2 data. This is because we in- 
clude both QCD evolution and higher twist in our 
fits, while QCD evolution is completely neglected 
in these other analyses || . We also compare the 
predictions with the standard GRV94 PDFs and 
with our modified GRV94 PDFs (LO+HT) to a 
few representative high energy CCFR and 77„ 
charged-current differential cross sections jn],[| 
on iron (these were not included in our fit). 
In this comparison we use the PDFs to obtain 
Fi and xF 3 and correct for nuclear effects on 
iron 1 1 |] . The parameterization for this nuclear 
effect is shown in appendix A of this article. The 
structure function 2xF± is obtained by using the 
Rworid fit from reference , as discussed in ap- 
pendix A. There is very good agreement with the 
neutrino data. 

According to Bloom-Gilman |l^] duality and if 
we use the x w scaling variable, the PDFs should 
also provide a reasonable description of the aver- 
age value of the structure functions in the reso- 
nance region. Figure ^ shows a comparison be- 
tween resonance data (from SLAC and Jeffer- 
son Lab, or parametrizations of these data [fl3|| ) 
and the predictions with the standard GRV94 
PDFs (LO) and with our modified GRV94 PDFs 
(LO+HT). There is good agreement with SLAC 
and JLab resonance data down to Q 2 = 0.07 (al- 
though these data were not included in our fit). 
In the Q 2 = limit, e.g. for E CM = 2 GeV, 
the modified PDFs yield a photoproduction cross 
section of 0.122 x 0.3/0.188 = 0.18 mb, which is 
in good agreement with experimental data. 

In order to have a full description of all charged 
current and processes, the contribution 
from quasielastic scattering jf4| must be added 
separately at x = 1. One may chose to use these 
modified GRV94 PDFs (LO+HT) to provide a 
description of all remaining inelastic scattering 
processes, including the resonance region, or one 
may chose to use them only above a certain value 
of invariant mass W, and add the lower mass res- 
onances separately. Probably the best prescrip- 
tion is to use these fits above the first resonance 
(e.g. above W=1.35 GeV) and add the contri- 
butions from quasielastic and first resonance |l5f 
(W=1.23 GeV) separately. In the resonance re- 
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Figure 4. Comparison of a SLAC and JLab low 
energy electron scattering data in the resonance 
region (or fits to these data) and the predictions 
of the GRV94 PDFs with (LO+HT, solid) and 
without (LO, dashed) our modifications. 



gion at higher mass (e.g. W=1.7 GeV) there is 
a significant contribution from the deep-inelastic 
continuum which is not well modeled by the ex- 
isting fits jl5| to neutrino resonance data (and 
using PDFs is better). Note that for nuclear tar- 
gets, nuclear corrections must also be applied. 

In conclusion, we present the result of a first at- 
tempt at modifying GRV94 PDFs (LO+HT) such 
that they provide a good description of e, [i and 
Up data both at low and high energies (including 
the average of the cross section in the resonance 
region). We have initiated a collaboration with 
scientists at Jefferson Lab to institute further im- 
provements such as allowing for different higher 



twist parameters for u, d, s, c quarks (e.g. hydro- 
gen versus deuterium and valence versus sea) . In 
addition, one can multiply the PDFs by a modu- 
lating function [§ A(W,Q 2 ) to improve modeling 
in the resonance region and comparing to Jeffer- 
son lab data [ fl3[ . In addition, we plan to include 
data on deuterium and heavier nuclear targets. 
Note that because of the effects of experimental 
resolution and Fermi motion (for nuclear targets), 
a description of the average cross section in the 
resonance region may be sufficient for some neu- 
trino experiments. It is expected that there are 
some differences between the form factors of reso- 
nances for and e/fi scattering. As a test of our 
approach, we will also compare our predictions for 
scattering to measured cross sections in the 
region of first resonance Jig ], where the largest 
differences are expected. 

A. Appendix 

In leading order QCD (e.g. GRV94 LO PDFs), 
F 2 for the scattering of electrons and muons 
on proton (or neutron) targets is given by the 
sum of quark and anti-quark distributions (each 
weighted the square of the quark charges) : 

F 2 (x) = X^e 2 [xqi(x) + xq^x)] . (I) 

Thus, our modified F 2 to describe low energy data 
down to photoproduction limit is given by 

O 2 

F 2 (x) = -^—^ F 2 (x w ), (2) 

V ; Q 2 + 0.188 1 h 1 ' 

where x w = x{Q 2 + 0.624)/(Q 2 + 1.735a;). 

In neutrino scattering we use the same modified 
scaling variable and the same correction factor in 
F 2 and xF 3 . 

In the extraction of original GRV94 LO PDFs, 
no separate longitudinal contribution was in- 
cluded. The quark distributions were directly fit 
to F 2 data. A full modeling of electron and muon 
cross section requires also a description of 2x F\. 
We use a non-zero longitudinal R in reconstruct- 
ing 2xF\ by using a fit of R to measured data. 
Thus, 2x Fx is given by 

2xF 1 {x) = F 2 {x)(l + 4Mx 2 /Q 2 )/(l + R). (3) 

In reality, a reconstruction of 2xF\ from the 
values of R and F 2 in neutrino scattering is not as 
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simple as in the case of charged lepton scattering 
(because of charm production). For charm pro- 
duction, the Bjorken scaling variable x no longer 
represents the fractional momentum carried by 
the struck quark in the infinite momentum frame 
due to the non-zero heavy mass of the charm 
quark (m c ~ 1.3 GeV). For charm production 
processes, the variable x is replaced by the slow 
rescaling variable £ = (1 + m 2 c /Q 2 )x. Therefore, 
the structure functions for the charm produc- 
tion (cp) and non-charm production (ncp) com- 
ponents are given by the following expressions. 



pcp+ncp^ = 



piicp 



2xF cp+ncp 



(X) = 



4M 2 x 2 /Q 2 



(4) 



R nc P(x) 
- AM 2 i 2 /Q 2 



pncp 



(X) 



*r p (o (5) 



1 + R nc P(£) 

In this model, the R wor id fit jl(| is used for R ncp 
and R cp , R WO rM is parameterized by: 

0.0635 



R 



world 



(x, 



log(Q 2 /0Mf {XlQ ) 
0.5747 0.3534 



Q 4 + 0.09' 



(6) 



where 6=1.+ 



12Q 2 



0.125 z 



Q 2 + 1.0 " 0.125 2 +a: 2 - 

The R wor id function provides a good descrip- 
tion of the world's data in the Q 2 > 0.5 and 
x > 0.05 region (where most of the R data are 
available). Fig || shows a comparison of all avail- 
able data on R and R wor id- Rworid is shown as 
the dotted line. Because of the effect of the charm 
quark final state mass in charm production in 
neutrino scattering, the effective R for neutrino 
scattering in the very low x region ( R e ff - solid 
line) is somewhat higher than R wor id (which is 
a fit to electron and muon scattering data). We 
use Eq. 4 and 5 to construct R e ff for neutrino 
scattering. However, this effect is not important 
at low neutrino energies (which are below charm 
production threshold). We include this effect be- 
cause we want our formalism to be also correct at 
high energies (e.g. in the CCFR/NuTeV region). 

Note that the R wor id function breaks down be- 
low Q 2 = 0.3. Therefore, we freeze the function 
at Q 2 = 0.35 and introduce the following function 
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Figure 5. Comparison of R data and R wor id- The 
dotted line represents R WO rid (for electron/muon 
scattering). The solid line shows the Reff val- 
ues for neutrino scattering, reconstructed using 
Eq. 4 and 5. The predictions with NLO resum- 
mation |l(J and NNLO calculations jllj are also 
shown for comparison. 



for R in the Q 2 < 0.35 region. The new func- 
tion provides a smooth transition from Q 2 = 0.35 
down to Q 2 = by forcing R to approach zero at 
Q 2 = as expected in the photoproduction limit 
(while keeping a 1/Q 2 behavior at large Q 2 and 
matching to R WO rid at Q 2 — 0.35). 



R(x,Q 2 



Rworid 

3.207 x 

^ Rworid 



Q 2 



Q 4 + i 
(x,Q 2 -- 



0.35). 



(7) 



In neutrino scattering the value of R is required to 
approach zero at Q 2 = only for the vector part 




Figure 6. The ratio of F2 data for 
heavy nuclear targets and deuterium as mea- 
sured in charged lepton scattering experi- 
ments(SLAC,NMC, E665). The band show the 
uncertainty of the parametrized curve from the 
statistical and systematic errors in the experi- 
mental data [ttSfl- 



of the interaction. However, the overall contribu- 
tion from R is expected to be small in this region. 
Therefore, it is reasonable to use Eq. g for R 111 
both the electron/muon and neutrino scattering 
cases. 

In the comparison with CCFR charged-current 
differential cross section on iron, a nuclear correc- 
tion for iron targets is applied. We use the fol- 
lowing parameterized function, f(x) (fit to exper- 
imental electron and muon scattering data for the 
ratio of iron to deuterium cross sections, shown in 
Fig |^) , to convert deuterium structure functions 
to (isoscalar) iron structure functions fl8[| ; 

f(x) = 1.096 - 0.364x - 0.278e- 21 ' 94:E 

+2.772Z 14 - 417 (8) 

For the ratio of deuterium cross sections to 
cross sections on free nucleons we use the follow- 
ing function obtained from a fit to SLAC data 
on the nuclear dependence of electron scattering 
cross sections f|. 

f(x) = (0.985 ± 0.0013) x (1 + 0.422a; - 2.745a; 2 



Figure 7. The total correction for nuclear ef- 
fects (binding and Fermi motion) in the deuteron, 
F^/F^ , as a function of x, extracted from fits 
to the nuclear dependence of SLAC F2 electron 
scattering data. 



+7.570a; 3 - 10.335a; 4 + 5.422a; 5 ). (9) 

This correction shown in Fig. 7 is only valid in the 
0.05 < x < 0.75 region. In neutrino scattering, 
we use the same nuclear correction factor for F2, 
xF 3 and 2a; F\. 

The d/u correction for the GRV94 LO PDFs is 
obtained from the NMC data for F 2 D /F 2 P . Here, 
Eq. ^| is used to remove nuclear binding effects 
in the NMC deuterium F2 data. The correction 
term, 8{d/u) is obtained by keeping the total va- 
lence and sea quarks the same. 

8(d/u)(x) = -0.0161 + 0.0549a; + 0.355a; 2 

-0.193a; 3 , (10) 

where the corrected d/u ratio is (d/u)' = (d/u) + 
5(d/u). Thus, the modified u and d valence dis- 
tributions are given by 

(11) 



1 



d' 



d v + u v S(d/u) 



(12) 



The same formalism is applied to the modified u 
and d sea distributions. Accidently, the modified 



7 



u and d sea distributions (based on NMC data) 
agree with the NUSEA data in the range of x be- 
tween 0.1 and 0.4. Thus, we find that any futher 
correction on sea quarks is not necessary. 
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